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ABSTRACT 

We present here new results of two-dimensional hydrodynamical simulations of the 
I ' eruptive events of the 1840s (the great) and the 1890s (the minor) eruptions suf- 

Q I fered by the massive star rj Car. The two bipolar nebulae commonly known as the 

Homunculus and the little Homunculus were formed from the interaction of these 
eruptive events with the underlying stellar wind. As in previous work (Gonzalez et 
al. 2004a, 2004b), we assume here an interacting, nonspherical multiple-phase wind 
scenario to explain the shape and the kinematics of both Homunculi, but adopt a 
more realistic parametrization of the phases of the wind. During the 1890s eruptive 
event, the outflow speed decreased for a short period of time. This fact suggests that 
the little Homunculus is formed when the eruption ends, from the impact of the post- 
outburst 7] Car wind (that follows the 1890s event) with the eruptive flow (rather than 
by the collision of the eruptive flow with the pre-outburst wind, as claimed in previous 
models; Gonzalez et al. 2004a, 2004b). Our simulations reproduce quite well the shape 
and the observed expansion speed of the large Homunculus. The little Homunculus 
(which is embedded within the large Homunculus) becomes Rayleigh- Taylor unstable 
and develop filamentary structures that resembles the spatial features observed in the 
polar caps. In addition, we find that the interior cavity between the two Homunculi 
is partially filled by material that is expelled during the decades following the great 
eruption. This result may be connected with the observed double-shell structure in the 
polar lobes of the tj Car nebula. Finally, as in previous work, we find the formation 
' of tenuous, equatorial, high-speed features that seem to be related to the observed 

equatorial skirt of r/ Car. 

Key words: stars: individual (rj Carinae) — stars: winds, outflows — hydrodynamics 
— shock waves 
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1 INTRODUCTION 

Located at a distance of 2.3 kpc, one of the most massive 
stars in our Galaxy ever discovered, rj Car is a well-known 
example of the evolved and unstable luminous blue variable 
(LBV) stars, characterized by sporadic, violent mass-loss 
eruptive events (e.g., Humphreys & Davidson 1994). The 
19th century spectrograms of ry Car (see also Walborn & 
Liller 1977; Humphreys, Davidson & Smith 1999 and refer- 
ences therein) provide evidence that the star underwent a 
giant eruption during ~20 yr around the 1840s, where a few 
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solar masses of gas (JSlO Mq) were expelled into the inter- 
stellar medium, and the total luminous output of ^lO'*^'^- 
^10^'^ erg was comparable to a supernova explosion (Smith 
et al. 2003a). It is still not understood how and why this out- 
burst occurred, but it may be connected with rj Car being a 
binary star system (e.g., Damineli 1996; Damineli, Conti, & 
Lopes 1997). From this event, a symmetric, bipolar nebula 
known as the "Homunculus" (that extends from -8 to -1-8 
arcsec along its major axis) was produced (Humphreys & 
Davidson 1994; Currie et al. 1996; Davidson & Humphreys 
1997; Smith & Gehrz 1998). During the 1890s, the historical 
light curve of rj Car also shows a much fainter eruptive event 
of shorter duration (~10 yr), but still an outburst in the 
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sense that the mass loss rate was enhanced compared with 
the normal underlying wind. This eruption eruption resulted 
in the formation of a smaller nebula embedded within the 
large Homunculus (with an angular size of about ±2 arcsec) 
called the "Little Homunculus" (LH) (Ishibashi ot al. 2003; 
Smith 2005). With a total mass of ~0.1 M©, the polar caps 
of the LH currently move at slower speeds (~ 250 km s~^) 
than the large Homunculus in the polar direction (~650 km 
s~^), but shares the same prolate geometry. The total kinetic 
energy released in the 1890s event is ~ 10'"'"'' erg (a factor 
of ~ 10^ smaller than the kinetic energy expelled during the 
great eruption). Apart from both Homunculi, observations 
of Tj Car with the Hubble Space Telescope (e.g., Morse et 
al. 1998) also show the presence of an equatorial skirt (or- 
thogonal to the axis of the Homunculus) that moves with 
velocities of 100-350 km s"^ (Smith & Gehrz 1998; David- 
son et al. 2001), although high speed features, with typical 
velocities of ~750 km s~^ (Smith et al. 2003b) or even larger 
(Weis 2005) have also been detected. These equatorial com- 
ponents may have been expelled from both the 1840s and 
the 1890s outbursts (e.g., Davidson & Humphreys 1997). 
The large Homunculus is the central portion of a larger neb- 
ula (the outer cjecta) around rj Car that extends up to a 
diameter of 60 arcsec (e.g., Weis 2001, 2005; Walborn 1976; 
Walborn et al. 1978). The outer ejecta are not symmetric 
and contain numerous filaments and knots that may have 
been produced by previous stellar wind mass loss. The sizes 
and morphology of such structures in the outer ejecta are 
manifold. Radial velocities for the knots reach up to 2000 
kms"'^ (Weis 2001, Weis 2005), but the average expansion 
velocity of the outer ejecta lies at lower values (around 750 
kms"'^), similar to the speeds found in the Homunculus. The 
outer cjecta also follow a bipolar pattern, like the large Ho- 
munculus, with approximately the same ax;is of symmetry. 

At optical wavelengths, the large Homunculus is mainly 
a reflection nebula. Its spectrum shows the presence of dust 
scattered emission (e.g., Hillier & Allen 1992) that allows 
one to see indirectly the shape of the r) Car wind (see Smith 
et al. 2003b) and also low-excitation intrinsic emission. In 
addition, near-infrared spectra obtained by Smith (2006) 
confirmed the existence of a double-shell structure at the 
edges of the polar lobes of the Homunculus (which was pre- 
viously inferred from thermal dust emission; see Smith et 
al. 2003a). A thin outer shell of intrinsic H2 emission (that 
traces the main scattering layer seen in visual images), and 
a thicker inner skin of [Fe II] (which partially fills the inte- 
rior of the lobes). On the other hand, the outer ejecta are an 
emission nebula. Smith & Morse (2004) present optical spec- 
tra of rj Car showing strong oxigen lines in some emission 
features of the outer ejecta. At high energies, Chandra has 
detected X-ray omission from the large-scale nebula around 
the star (and also from the central object, probably due to a 
wind collision region around the massive binary system; e.g., 
Corcoran ct al. 2001; Pittard & Corcoran 2002). Soft X-ray 
(0.1-0.8 keV) observations by Seward et al. (2001) show an 
extended shell clearly associated with the debris field of the 
outer ejecta which are compatible with a collision process 
between a fast wind and cometary knots of slower material. 
The ri Car nebula is also a bright radio source consistent 
with thermal (free-free) emission. Observations by Retallack 
(1983) at 1.415 GHz (taken with the Fleurs synthesis tele- 
scope) show emission from an overall region as large as 40 



arcsec (i.e., larger than the optical Homunculus). Gonzalez 
et al. (2006) estimated the contribution by shocks to the to- 
tal radio-continuum emission detected from the Car neb- 
ula. Using observational estimates of the wind parameters 
of the eruptive event of the 1890s and of the stellar wind 
after the end of the eruption, these authors investigated the 
evolution of the polar caps of the LH formed as a result of 
the collision between these outflows. They found that the 
LH emits continuum radiation which is detectable at radio 
wavelengths and indeed, has an appreciable contribution to 
the total fiux of the r] Car nebula. 

Different models have been proposed to explain the 
shaping and kinematics of the rj Car bipolar nebulae (see, 
for instance, Soker 2001; Soker 2004; Matt & Balick 2004; 
Gonzalez et al. 2004a, b; A. Frank et al. 1995, 1998). One 
possible explanation is that it is produced by the interaction 
of the winds expelled by the central star at different injec- 
tion velocities (e.g., Icke 1988; Frank, Balick & Davidson 
1995; Dwarkadas & Balick 1998; Frank, Ryu & Davidson 
1998; Langcr, Garcia-Segura & Mac Low 1999; Gonzalez et 
al. 2004a,b). Adopting a colliding wind scenario, Gonzalez 
et al. (2004a,b) performed two-dimensional numerical simu- 
lations of the Homunculi nebulae of rj Car. In their models, 
the large Homunculus is formed by the interaction of the 
eruptive outflow of the 1840s with the pre-outburst r] Car 
wind (both with different degrees of nonspherical symme- 
try). A second eruption (assumed to be spherical) collides 
with the pre-outburst wind giving rise to the LH. These au- 
thors showed that such a scenario could explain the shape 
and kinematics of the Homunculi, and also the existence of 
the high-velocity features observed in the equatorial plane. 
However, their prescribed conditions to create the LH nebula 
do not agree with overall results suggest by the observations. 
During the 1890s event, the rj Car wind slowed down and in- 
creased its mass loss rate, returning to its normal quiescent 
state after the eruption ended. Another observational re- 
sult (as mentioned earlier) is that the Homunculus contains 
much more mass than that that had been previously recog- 
nized (Smith et al. 2003a). In this work, we perform new 
numerical simulations that incorporate these observational 
results and investigate the overall evolution of the bipolar 
outflows of T] Car. 

The paper is organized as follows. In § 2, we describe 
the model. The numerical simulations and the discussion of 
the results are presented in § 3, and in § 4 we draw our 
conclusions. 



2 THE MODEL 

In our model, we consider (see also Gonzalez et al. 2004a,b) 
a simplified interacting wind scenario, in which a nonspheri- 
cal outburst collides with a slow wind also with asymmetric 
density and velocity distributions. We suppose that rj Car 
originated a nonspherical environment (prior to the major 
eruption of the 1840s) from the ejection of a pre-outburst 
wind into a homogeneous ambient medium. According to 
observational estimates, we also assume that both the mass 
loss rate and the injection velocity were drastically increased 
during the main eruption of the 1840s, after which the orig- 
inal slow wind resumes. We also suppose that the wind 
parameters are suddenly changed during the smaller mass- 
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Figure 1. The solid line gives the best fitting curve to the mea- 
sured velocity in H2 (squares) by Smith (2006) at different lat- 
itudes of the 7] Car Homunculus (with 9= 0° at the pole). For 
this fitting we have assumed A = 1.9, di= 670 km s~^, V2 = 100 
km s~^, which gives (see eq. [T]) a polar expansion velocity Vp = 
657.53 km s~^ and Ve= 112.47 km s^^ at equator. 



ejection event of the 1890s, after which the original wind 
again resumes. However, observations of 77 Car at this epoch 
(e.g. Whitney 1952; Walborn & Liller 1977) give evidence 
that the expansion speed of the stellar wind decreased dur- 
ing this eruption, rather than increasing (see also Le Sueur 
1870). Such variabilities in the wind parameters of Car 
during the 19th century outbursts must have resulted in the 
formation of a pair of double-shock wave structures (called 
working surfaces; see Raga et al. 1990), that correspond - 
in our colliding wind scenario - to the large and little Ho- 
munculi. Then, the large Homunculus would be produced by 
the interaction between an outburst wind (with a latitude- 
dependent velocity and density) and a pre-eruptive slower 
wind (also with a nonspherical symmetry) . Nevertheless, the 
little Homunculus would be formed - unlike previous numer- 
ical modeling (Gonzalez et al. 2004a, b) - from the impact of 
the post-outburst wind with the outflow expelled during the 
minor event, that is, when the eruption ended. According to 
Gonzalez et al. (2006) , a sudden increase in the ejection ve- 
locity at the end of the eruption instantaneously occurs and 
forms the inner nebula at the base of the wind (as the fast 
upstream flow begins to be ejected). 

In order to estimate the flow parameters (such as the in- 
jection velocity v and the number density n) for the different 
wind phcises, we have fitted a curve to the expansion speed 
of the outer H2 shell of the rj Car Homunculus measured by 
Smith (2006). From this fitting, we find a latitude-dependent 
velocity given by. 



Figure 2. Behaviour of the adopted parameters of the interact- 
ing outfiows. We show the mass loss rate (top panel) and the wind 
velocity at the poles (bottom panel) as functions of time. During 
the great eruption of the 1840s (~20 yr), both parameters were 
suddenly increased, while during the minor event of the 1890s 
the flow velocity decreased for ~10 yr. In addition, we assume a 
gradual fading in the mass loss rate during the years after the end 
of both events, asymptotically approaching the mass loss rate of 
the underlying wind (10""^ Mq yr~^). 



with 



F{0) = 



{V2/V1) + 



1 e2 ^ 

where the parameter z = Acos(20) controls the shape of the 
Homunculus, A is a constant and 9 is the polar angle. The 
velocities vi and V2 are related to the expansion speed in the 
polar {6= 0°) and equatorial {9= 90°) directions. In Figure 
1, we present the best fit to Smith (2006)'s observations 
obtained with A = 1.9, 111= 670 km s~^, «2= 100 km s~^. 
With these parameters, the predicted expansion velocity in 
the polar and equatorial directions are Vp— 657.53 km s~^ 
and Ve— 112.47 km s~^, respectively. Using this fit of the 
expansion velocity of the outer Homunculus, we estimate 
the flow parameters of the interacting winds as follows. 

From a straightforward application of the formalism 
developed by Canto et al. (2000) for outflows with time- 
dependent injection velocities (see also Gonzalez & Canto 
2002; Canto et al. 2005), we studied the dynamical evolu- 
tion of the Homunculus. Let us first consider the downstream 
wind with latitude-dependent injection velocity vq (given by 
eq. [1]) and constant mass loss rate per unit solid angle. 
Then, the injection number density of the wind is equal to, 



(2) 



v = vi F{9), 



(1) 



where no{— rno/4:TVfivorQ; being fi the mean mass per nu- 
cleon and rho the total mass loss rate) is the number density 
at the injection radius ro (assumed to be a few stellar radii). 
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TABLE 1 



Parameters of the colliding outflows 



Wind Phase 


A 


vi [km s -"-] 


1)2 [km s -"-] 


Vp [km s ■"-] 
(e = 0°) 


tjg [km s ■"-] 
(6» = 90°) 


m [Mq yr-i] 


Pre- outburst wind 


2.4 


250.00 


14.00 


248.07 


15.93 


10-3 


Great Eruption 


1.9 


687.76 


102.65 


674.95 


115.45 


5 X 10-1 


Post-outburst wind 


1.9 


500.00 


14.00 


489.37 


24.63 


10-3 X (t) 


Minor Eruption 


1.9 


200.00 


10.00 


195.84 


14.16 


10-2 


Post-outburst wind 


1.9 


500.00 


300.00 


495.62 


304.38 


10-3 X ^{t) (t) 



(t) see eq. [6] 
(t) see eq. [7] 



Note that F{9) depends on the particular values of «i, V2, 
and z of the pre-outburst wind. 

When the great eruption begins, the wind parameters 
are suddenly increased to avo and brho, respectively (where a 
and b are constants). Such a variation in the flow parameters 
forms instantaneously (at the base of the wind) a working 
surface (which will correspond in our model to the large 
Homunculus) that moves with a constant velocity given by, 



(3) 



with a = {a^'^ + ah^^"^) / {a^'^ + 6^/^). This velocity is in- 
termediate between the low- velocity downstream wind (wo) 
and the faster upstream outflow {avo). 

Given (for instance) vo and v-ws , equation (|3} gives a as 
a function of a and 6. We estimate b from observations of the 
mass contained in the Homunculus and the duration of the 
major eruption (and compare with the mass loss rate of the 
standard wind of ri Car) and then, we obtain a. Assuming 
the same latitude-dependent velocity and density (eqs. [T] 
and [2] , respectively) for the pre-outburst wind and the great 
eruption, and using the fit to Smith (2006) 's observations of 
the expansion velocity of the Car Homunculus (Fig. 1), we 
finally obtain the flow parameters at the different phases of 
the interacting winds, keeping vq as a free parameter. 

The historical light curve of Car shows a sudden in- 
crease in brightness when the eruptions of the 1840s and 
1890s turned-on and a gradual fading during the years after 
the end of the eruptive events (e.g. Walborn & Liller 1977; 
Davidson 1987; Humphreys & Davidson 1994; Humphreys, 
Davidson & Smith 1999). This suggests that the post- 
outburst mass-loss rates of the major (mi) and the minor 
(7712) events decrease, approaching the mass loss rate of the 
underlying wind (mo = 10-3 M0 yr-^; e.g. Humphreys & 
Davidson 1994; Davidson & Humphreys 1997). We then con- 
sider, 



mi = mo (t){t), 
and 

m2 = mo <p{t), 
respectively, with 



(4) 



(5) 



and. 



+ 1 



mo 



mo 



TV I t 
2 



ti 



ifi{t) = — h I 1 : — I sm 



mo 



mo 



TT I t 
2 



At2 



(6) 



(7) 



where rhge{— 5x IQ-'^ Mq yT~ ) and rhme{— 10~ Mq yT~ ) 
are the estimated mass loss rates during the eruptions, ti 
and t2 correspond to the transition times when both events 
end, and Ati{= 30 yr) and At2{= 10 yr) represent in our 
model the duration of the post-eruption phases. In Figure 
2, we show the behaviour of the adopted parameters (the 
mass loss rate and the ejection velocity at the poles) of the 
interacting outflows as functions of time. 



3 NUMERICAL SIMULATIONS 

We have perfomed gasdynamic 2D numerical simulations 
(considering axial symmetry) of the 19th century outbursts 
of rj Car. Adopting the colliding wind model described in 
§2, we use the adaptive-grid YGUAZU-A code originally 
developed by Raga et al. (2000; see also Raga et al. 2002) 
and modified by Gonzalez et al. (2004a, b). This code in- 
tegrates the hydrodynamic equations explicitly accounting 
for the radiative cooling with a set of continuity equations 
for the atomic/ionic species HI, HII, Hel, Hell, Helll, CII, 
CHI, CIV, 01, OH, and OIII. The fiux-vector splitting algo- 
rithm of Van Leer (1982) is employed. The simulations were 
computed on a five-level binary adaptive grid with a maxi- 
mum resolution of 3.9 x 10^* cm, corresponding to 1024 x 
1024 grid points extending over a computational domain of 
(4 X 10^^ cm) X (4 X 10^^^ cm) . The adopted abundances (by 
number) for the different elements are (H, He, C, O) — (0.9, 
0.099, 0.0003, 0.0007). 



3.1 Initial Physical conditions 

In Table 1 (see also Fig. 2), we list the adopted parame- 
ters for the interacting winds. We have assumed different 
equator-to-pole velocity (and density) contrasts. 
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Figure 3. Stratifications of temperature (top left), density (top riglit), pressure (bottom left), and velocity-field superposed on the 
density map (bottom riglit) obtained 169 yr after the great eruption of 7] Car (i.e., around 2009) are presented. The arrow shown at the 
lower right corner of the velocity-field stratification corresponds to 1000 km s~^ (see the text for more details) 



Initially, the computational domain is filled by a homo- 
geneous ambient medium with temperature Ta = 10^ K and 
density Ua = 10"'^ cm""*. Afterwards, a gaseous toroidal dis- 
tribution is formed by the injection of a pre-outburst wind 
into this unperturbed environment (with a different A pa- 
rameter from the other injected later wind phases which 
controls the shape of the wind). For this wind, we have 
assumed a terminal velocity (in the polar direction) vi = 
248.07 km (Gonzalez et al. 2004a, b) and a mass loss 
rate of 10~^ M© yr~^ which was injected at a distance of 
ro=10^®cm (from the stellar surface) with a temperature of 
ro=10^ K. Once a poloidal environment is produced, a much 
faster {— 674.95 km s~^) and massive outflow (5 xlO~^ Mq 
yr~^) is expelled during the estimated duration of the great 



eruption (~ 20 yr; e.g. Davidson & Humphreys 1997). For 
this outburst phase (that dominates the momentum flux) 
we have adopted the same degree of nonspherical symmetry 
as the observed for the large Homunculus, as given in Figure 
1. After this event, a third outflow (the post-outburst wind) 
is turned on with similar conditions to the current stellar 
wind of T] Car (e.g. Le Sueur 1870; Smith et al. 2003b), that 
is, we assume an injection velocity of 489.37 km s^^ at the 
poles and a time-dependent mass loss rate that approaches 
10""^ M0 yr~^ (eqs. [J and [B]). In order to account for the 
minor eruption, during 10 years (and 50 years after the great 
eruption event), we assumed a slower eruptive wind with a 
terminal velocity - along the symmetry axis - of ~ 195.84 
km s~^ and a mass loss rate of 10~^ M© yr~^. After this 
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outflow, a faster 



495.62 km 



but less massive wind 



with a mass loss rate approaching to 10~^ M0 yr~^ (eqs. 
[5] and [7]) resumes. We also have assumed for this wind a 
higher speed {= 304.38 km s~^) at equator, which is con- 
sistent with the present-day latitudinal structure in r] Car's 
stellar wind (Smith et al. 2003b). 

We should notice that Smith (2006) calculated the mass 
loss distribution assuming uniform density as a function of 
latitude with a constant width. These are not the assump- 
tions of our current models. Instead, we assumed that both 
the ejection velocity and the injection density depend on 
latitude (eqs. [1] and [2], respectively) and therefore, a con- 
stant mass loss rate per unit solid angle is adopted (see §2). 
However, a similar scenario to the Smith's findings was pre- 
viously study in Gonzalez et al. (2004b), where we assumed 
a nonspherical outburst wind of the 1840s impinging on a 
slow pre-outburst wind with a larger mass-loss rate in the 
polar direction (run D). We found that this scenario does 
not develop significant equatorial features. 



g 
3, 



3.2 Results of the simulations 

In this section we present the two-dimensional hydrodynam- 
ical numerical simulations performed for the outbursts of rj 
Car, adopting the colliding wind scenario described in § 2. 
Figure 3 shows the temperature, density, pressure, and ve- 
locity maps computed for the interaction of the five winds 
above at a time t= 169 yr of evolution after the great erup- 
tion. As predicted in §2, the simulations show the formation 
of the outer Homunculus with a double-shock structure, hav- 
ing an inner shock that decelerates the fast outburst fiow and 
an outer shock that accelerates the lower velocity precursor 
wind. At the poles, a cold thin shell behind the inward shock, 
and a hotter and thicker region behind the outward shock 
are formed. This difi'erence is due to the radiative cooling of 
the shocked material behind both shocks (see Gonzalez et 
al. 2004a). We find that for the great eruption, the momen- 
tum flux is dominated by the material expelled during the 
eruptive event (that is, the pre-outburst wind to the great 



eruption momentum ratio, mo vo / nig 



<^ 1 (where rhg 



and Vge are the mass loss rate and the wind velocity during 
the eruption, respectively), so that the pre-outburst wind 
has a minor effect on the kinematics and morphology of 
the bipolar lobes. Since the mass loss rate during the erup- 
tion is more than two orders of magnitude (a factor of 500; 
see Table 1) larger than that of the pre-outburst wind, the 
Homunculus retains almost the same degree of asymmetry 
imprinted in the 1840s event near the star. 

We note that the shape and kinematics of the outer ex- 
panding shell resembles that of the large Homunculus previ- 
ously simulated by Gonzalez et al. (2004a, b) , but important 
differences from these previous models are identified in the 
embedded structures of Figure 3. In the case of the minor 
eruptive event, the momentum flux is also dominated by the 
eruption, but the post-outburst wind has a signiflcant effect 
on the kinematics and morphology of the inner Homunculus 
once the low-velocity downstream material is completely in- 
corporated to the layer. Gonzalez et al. (2006) have shown 
that this happens at a critical time tc = (Ai) /(cr — 1), being 
At (=10 yr) the duration of the eruptive event. Adopting the 
flow parameters given in § 3.1 (cr ~ 1.25; a= 2.53 and b=0.1), 
we obtain tc — 39.4 yr. In FigureH] we show the behaviour of 



t LyrJ 

Figure 4. Behaviour of the velocity of the polar caps of the 
inner Homunculus as a function of time. Initially they move with 
constant velocity (~ 245.55 km s~^). At a critical time t= 39.4 
yr (when the downstream shock disappears), the caps begin to be 
accelerated, asymptotically approaching the velocity of the post- 
outburst wind (~ 495.62 km s~^). After 109 yr of evolution (i.e., 
around the year of 2009) they expand at ~ 297 km s"!. This 
value is consistent with the observed current expansion speed of 
~ 250 km s-i of the little Homunculus (Smith 2005). 



the velocity of the polar caps of the inner Homunculus as a 
function of time. Initially, they move with constant velocity 
(~ 245.55 km s~^; see eq. [3]), until the downstream shock 
disappears and a one-shock structure stage begins. Later, 
the polar caps begin to be accelerated, asymptotically ap- 
proaching the velocity of the post-outburst wind (see also 
Gonzalez & Canto 2002). After 109 yr of evolution (around 
year 2009), they expand at ~ 297 km s~^ and are located at 
a position Vs — 9.07x lO'^® cm from the star. At a distance 
of 2.3 kpc, rs corresponds to an angular size of ± 2.6 arcsec, 
which is consistent with the angular extent of roughly ±2 
arcsec (along the major axis) of the inner Homunculus mea- 
sured by Ishibashi et al. (2003). This is also consistent with 
the angular size (± 3 arcsec) observed in the simulations. 

The inner Homunculus quickly becomes Rayleigh- 
Taylor unstable (see flgure[5]) due to the interaction of the 
low-density fast wind (post-outburst wind) that pushes and 
accelerates the high-density slow wind (minor eruption) . The 
growth time for this instability can be estimated as follows. 

The density of the fast wind at the position of the shock, 
Vs, is obtained by. 



Pfw — 



(8) 



where Mf^ is the mass loss rate of the wind, and Vf^ is its 
velocity at r^. On the other hand, the shock-bounded layer 
of the slow wind has a surface density asm given by. 



(9) 
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Figure 5. Detail of tlie density map of Fig. 3, showing the lit- 
tle Homunculus at t= 119 yr after the minor eruptive event of 
Car (or around the year 2009). The low-density fast wind (post- 
outburst wind) pushing the high-density slow wind (minor erup- 
tion) causes the polar caps of the little Homunculus to become 
Rayleigh- Taylor unstable. 

where Maw = Maw At is the mass in the shock-bounded 
layer, At is the time interval of injection of the slow wind, 
and Maw is its mass loss rate. In the above, we have assumed 
that the shock-bounded layer is thin. 

Initially, the shock-bounded layer (and the slow wind 
shock) moves with a velocity Va given by (eq. [3] ) , 

Va — UVaw (10) 

where Vaw is the velocity of the slow wind. In the shock 
reference frame, the fast wind moves with velocity v = Vfw — 
Va, and thus excerts a hydrodynamical pressure pfwV^ on the 
shock-bounded layer. This layer will therefore, experience an 
acceleration as — PfwV^ /o'aw- Taking the values of the mass 
loss rates and wind velocities from Table 1, we obtain a 
value for the acceleration experienced by the layer of a^ = 
3.96 xlO"^ cm s"^. 

The dispersion relation for the Rayleigh- Taylor (RT) 
instability is calculated by. 



p2 + pi 

where p2 > pi are the densities of the fluid layers, g is the 
acceleration, and is a wave number. We identify p2 = paw , 
Pi = Pfw and g — aa. Considering that the growth rate for 
the RT instability is largest for the shortest wavelengths, 
we adopt the minimum resolved wavelength as 3 times the 
resolution of our simulation, \rt = 1.17 x 10^^ cm. With 
Pfw = 10"^** g cm"^ and paw = 2 x 10"^'' g cm"^, the RT 
instability grows in about one year, in agreement with the 
simulations. 

In Figure 6, we show that our model of interacting non- 
spherical winds can explain the formation of the large and 
little Homunculi. From the simulations, we have found an 




Figure 6. Present-day structure of the rj Car nebulae. The in- 
ner and outer Homunculus and a tenuous equatorial skirt are 
depicted. A thick dense layer between both Homunculi is also 
observed. 



appropriate combination of the flow parameters (which con- 
trol the degree of asymmetry) that best matches the kine- 
matic structure of the thin outer shell of H2 observed by 
Smith (2006). In contrast with our previous models of ry 
Car (Gonzalez et al. 2004a, b) in which we observed a hollow 
nebula, it is noteworthy that the new numerical experiments 
show a thicker layer (with a density between ~ 10~^° g cm~'^ 
and 10~^^ g cm~^ along the symmetry axis) partially filling 
the interior of the lobes of the large Homunculus. It may be 
related with the inner skin of [Fe II] emission detected by 
Smith (2006). 

On the other hand, as in the previous work, the high- 
velocity ejections at equator, which arise from the impact 
of the external Homunculus on the shock front of the pre- 
emptive wind, resemble the observed equatorial skirt of r; 
Car. At a time t ~ 100 yr after the great eruption turn-on, 
the impact at low latitudes occurs producing a hot (T ~ 10^ 
K), tenuous (p ~ 10~^* g cm~^) structure. This structure 
is accelerated, reaching an expansion velocity of ~1000 km 
s~^ and produces fast ejected material that could be related 
to the observed material confined to the equatorial plane 
of 77 Car. The gas is concentrated to the equator by the 
impact of the two shock fronts, and then it flows into an 
environment of decaying pressure. Why are these equato- 
rial features moving faster and what accelerate them in our 
simulations are still open issues which we will attend in a 
forthcoming paper. This effect has been also addressed in 
Gonzalez et al. (20Q4a,b) and it is probably connected with 
a nozzle problem. As can be seen in the simulations (Figs. 
3 and 6), the equatorial ejection has some resemblance with 
a supersonic flow through a tube of varying cross-sectional 
area (the problem of the de Laval nozzle). At supersonic 
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speeds, an increase in velocity is expected when the area of 
the nozzle increases, such as our models show. 



4 DISCUSSION AND CONCLUSIONS 

In this article, we carried out high-resolution two- 
dimensional gasdynamic simulations of the dramatic out- 
bursts suffered by the star r] Car during the 1800s, the larger 
of which occurred in the 1840s and resulted in the forma- 
tion of the large Homunculus nebula and the smaller one in 
the 1890s created the inner little Homunculus. During these 
events, the parameters of the r) Car wind may have dras- 
tically changed in short periods of time. In contrast with 
the great eruption (where both the mass loss rate and the 
ejection velocity were suddenly increased), during the mi- 
nor event of the 1890s (much fainter and of shorter duration 
than the 1840s event) the mass loss rate was enhanced, but 
the wind velocity diminished compared with the normal rj 
Car wind. 

Considering a simplified interacting stellar wind sce- 
nario, we could explain the shape and the observed expan- 
sion speed (as a function of latitude) of the large Homuncu- 
lus. In addition, our numerical models show that the little 
Homunculus is formed at the end of the 1890s eruption, 
when the post-outburst 77 Car wind collides with the erup- 
tive outflow. Important differences with regard to our pre- 
vious models of the 77 Car nebula (Gonzalez et al. 2004a,b) 
have been obtained. In Gonzalez et al. (2004a,b), we as- 
sumed that the impact between the 1890 outburst and the 
pre-outburst wind was the cause of the development of the 
little Homunculus nebula. As a matter of fact, although the 
momentum flux is dominated by the 1890s eruption, the 
post-outburst wind has an important effect on the kinemat- 
ics and morphology of the little Homunculus as it causes 
the complete deposit of the material expelled during the 
eruption into the polar caps. At this time, the powerful stel- 
lar wind accelerates the inner Homunculus material so that 
it asymptotically reaches the velocity of the post-outburst 
wind. Due to Rayleigh- Taylor instabilities generated by the 
low-density fast wind pushing the high-density slow wind, 
the polar caps of the inner Homunculus quickly develop fil- 
amentary structuring (Fig. 5) that shows some resemblance 
with the observed spatial structures in the polar lobes of the 
little Homunculus by Smith (2005) that suggest that it is not 
perfectly homologous to the larger Homunculus nebula. 

Also, in contrast with the previous models of r] Car 
(Gonzalez et al. 2004a, b) , it is noteworthy that in the present 
numerical modeling the interior cavity between the Homun- 
culi is partially filled by material that is expelled during 
the decades following the end of the 1840s great eruption, 
rather than being almost empty. This agrees with the ob- 
served double-shell structure observed in the polar lobes by 



We notice that the effect described her is distinct from the one 
discussed, e.g., in Akashi & Soker (2008). These authors claim 
the presence of jets (or coUimated fast winds) launched by the by 
the central star (or by a companion) to explain the presence of 
an expanding disk and then relate its acceleration to compression 
of the material in the equatorial plane by the two lobes. 



Smith (2006), consisting of a thin outer H2 skin (which con- 
tains most of the material of the nebula) and a thicker [Fe 
II] layer with a more irregular spatial distribution. 

As in the previous models (Gonzalez et al. 2004a,b) , the 
present results show the formation of an equatorial outflow 
with both low and high velocity features. These are probably 
related to the equatorial skirt of r) Car. We note however, 
although it has been predicted that the observed skirt might 
be associated with the two outbursts (e.g. Davidson et al. 
2001), our simulations indicate that only the great eruption 
contributes to its formation. 

A final remark is in order. One could ask whether one 
can learn about the progenitor of Eta Car from the assumed 
mass loss history or what can cause such sharp changes in 
the mass loss rate and speed, without affecting the gen- 
eral shape of the Homunculi. As stressed in Gonzalez et al. 
(2004a,b). these questions are vey probably connected with 
Eta Car being a binary star system and the nature of the 
interaction between the main and companion stars. In this 
work, however, we have focussed in the formation and dy- 
namical evolution of the shock structures associated to the 
large and little Homunculi given the parameters (terminal 
velocity, mass-loss rate and/or density) of the different wind 
phases, without addressing the inner mechanism that flrst 
triggered it or its variability at the much smaller scales of the 
central source of Eta Car. The orbit of the secondary star, 
for instance, has apoastro and periastro distances of ~14 
AU and ~3 AU, respectively, while in our simulations, the 
outflows are injected at much greater distances (10^® cm). 
A detail study univocally connecting both the inner source 
scales and the outer scales was, therefore, out of the scope 
of the present study. Nonetheless, there have been some flrst 
efforts in this direction (e.g., Soker 2001 ; Falceta-Gongalvez 
& Abraham 2009) that should be further explored in the fu- 
ture. 
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